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Abstract 

The determination of smuon and neutralino masses in smuon pair production is an important 
part of the program of spectroscopic studies of Supersymmetry at a high energy linear collider. 
In this note we report the first results of a study of e^e^ — ^ P-rP-r ^ high-mass, cosmology- 
motivated Supersymmetric scenario at 3 TeV at CLIC. This process is a good example to study 
requirements on the beam energy spectrum and polarisation and the track momentum resolution 
in a simple final state. We discuss the expected accuracy on the mass measurements as a func- 
tion of the momentum resolution, luminosity spectrum, beam polarisation and time stamping 
capability. Results obtained at generator level ai^e validated by comparison to full simulation 
and reconstruction. Preliminary requirements for the detector performances and beam polarisa- 
tion are presented. 



1 Introduction 



One of the main objectives of linear collider experiments is the precision spectroscopy of new particles 
predicted in theories of physics beyond the Standard Model (SM), such as Supersymmetry (SUSY). Since 
some, or most, of these particles may have masses of ^(1 TeV), these studies may be central to the physics 
program of a multi-TeV e^e^ linear collider, such as CLIC. 

In this note, we study the production of the supersymmetric partners of the muon in a specific SUSY 
scenario, where we assume R-parity conservation within the so-called constrained Minimal Supersymmetric 
extension of the SM (cMSSM). In this model the neutrahno (^f) is the lightest supersymmetric particle and 
the specific parameters of the benchmark point |T| are chosen to make it compatible with current collider and 
cosmology data. In particular, the properties of the lightest neutralino are such that it generates the correct 
amount of relic dark matter density in the universe, as obtained from the analysis of the WMAP data lH. 
Scalar muons {jX^ and fl^) are the supersymmetric partners of the right- and left-handed charged muons. 
Smuons are produced in pair through ^-channel 7/Z exchange in the process e^e^ — > jUr jUr ^^id each decay 
into an ordinary muon and a neutralino, X\- The neutralino, being weakly-interacting, escapes detection. 
Therefore, the experimental signature of the process is two oppositely charged muons plus missing energy. 
This study concentrates on the lightest smuon, jX^, which, for the chosen model pai^ameters, has a mass of 
1108.8 GeV , while the mass of the lightest neutralino is 554.3 GeV. The accurate determination of their 
masses is an essential part of the spectroscopy study of a high mass SUSY scenario at CLIC. A study of the 
variation of the predicted relic dark mass density in the universe Q.h^ with the lightest neutralino mass in the 
cMSSM shows that a ±1.0 GeV uncertainty on its mass con^esponds to a ±0.05 relative uncertainty on Q.h^, 
i.e. the current accuracy from cosmic microwave background observations [3]. The main aim of this study 
is to assess the requirements for a detector at CLIC operating at a centre-of-mass energy, ^/s, of 3 TeV as a 
function of the track momentum resolution, luminosity spectrum and beam polarisation. The reconstruction 
of the particle masses thr^ough the endpoints of the muon momentum spectrum is a good example for these 
requirements since the analysis is particularly simple and can be canied out using a simple momentum 
smearing on generator-level observables. Results are validated using full simulation and reconstruction with 
the CLIC-ILD detector model. 



2 Simulation data sample 

The simulation is performed for the cMSSM parameters of point K' of ref. fl]. In the cMSSM the mass pa- 
rameters are defined at the GUT scale. The subsequent evolution to the electro-weak scale is performed 
using the renormalisation group equations of ISASUGRA 7.69 El. Signal events ai^e generated using 
Pythia 6.125 lHJ. At 3 TeV, the production cross section for the process e^e^ — )■ fi^fin, for unpolarised 
beams, is 0.71 fb. Beamstrahlung effects on the luminosity spectrum are included using results of the CLIC 
beam simulation for the 2008 accelerator parameters 10. Initial state radiation (ISR) is included in the 
event generation in Pythia. The following background processes have been included in the background 
calculation: 

Process Cross section 

e+e- W+W- ^ jx+ix-v^v^ a=10.5 fb 

e+e- tPz^ li+ii-vv a=0.5 fb 

e+e- Inclusive SUSY ^n+^-X a=0.4 fb 
e+e^ — > jUVg/iVg (inclusive SM) <7=135 fb 

The first three processes have been simulated with PYTHIA. In addition, the inclusive SM process 
e+e^ — > ix+lx^VeVe is generated using CompHep Q, removing the contributions from the VK+W^ and 
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Z*^Z^ diagrams, to avoid double counting. The estimated cross section is 135 fb. In the background study 
we neglect the e^e^ — /x+/i^v^v^ contribution not due to W+W^ and Z°Z° decays, since its cross section 
is only ~0.2 fb. We assume a data sample corresponding to an integrated luminosity of 2 ab^^ taken at 
a nominal ^/s energy of 3 TeV, corresponding to ~3.5 years (1 year = 10^ s) of run at the nominal CLIC 
luminosity of 6x10^^ cm^^s^^ Beam polarisation is in general extremely helpful in the study of SUSY 
processes both to improve the signal-to-background ratio and as an analyser |J8|. We consider here three op- 
tions for beam polarisation: i) unpolarised beams, ii) P(e^)=-i-80 % and P(e+)=0 % and iii) P(e^) = -1-80 % 
and P(e+)=-60 %. The main benefits of beam polarisation for this analysis ai^e the suppression of the W^W^ 
background (by a factor of five for ii, and ten for iii)), the enhancement of the smuon cross section (by a 
factor 1.5 for ii), and 2.3 for iii)) and the possibility to disentangle pLRpLR from pLujii and IXlP-l production. 
In this analysis we use observables at the generator level applying a track momentum smearing. Results are 
validated by comparing with fully simulated and reconstructed events in section 3.3.7. 

3 Analysis Procedure 

3.1 Signal topologies and event pre-selection 

The signal process has two undetected ^{''s in the final state. Therefore, the main characteristics of signal 
events are large missing transverse momentum, missing energy and acoplanarity (see Figure[T]). Despite the 




Figure 1: Display of a simulated e^e l^Ri^R l^-^l^ XiM 

striking signature of two muons and large missing energy, the small anticipated signal production cross sec- 
tion at the K' benchmark point, makes this analysis rather challenging. In our analysis, the signal selection 
proceeds as follows. First, we apply an event pre-selection, which requires two oppositely-charged muons 
with pt > 5 GeV and |cos0| < 0.985, where 6 is the particle polar angle. Next, we combine the values of 
the signal probabilities for the following discriminating variables into a global likelihood variable Prob: 

• visible energy E^is, 

• missing transverse energy Ej^,„iss, 

• sum of transverse momentum of the muons L |Pf I > 
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• maximum acoUinearity and acoplanarity, 

• polar angle of the missing energy vector (dmiss) 

• invariant mass of the two muons, 

• the thrust of the two muons, 

• unbalance of the muon momenta A 

• missing mass Mmis 

where A = ^1 — jp^^T^p^jr) andM^is = {s+M^^^ — ly/sE^isY^^ with the missing mass calculated from the 

visible energy f'vis and momentum /Vis, andMyis = (E^i^ ~^vis)'^^- Figures |2] and [3] show the distributions of 
some of the discriminating observables for signal and background samples after pre-selection and requiring 
> 2500 GeV. 
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Figure 2: Discriminating variables used in the combined likelihood function: (upper left) Eyis visible energy, 
(upper right) ^tmiss missing transverse energy, (lower left) ^pt sum of the pt of the muons and 
(lower right) invariant mass of the two muons 



3.2 Final selection efficiency and background estimate 

The normalised signal-to-background ratio, S/B, values of these variables, as well as the combined probabil- 
ity Prob are computed for different detector resolution assumptions: 5pt/p^= 2 x 10^^, 4 x 10^^, 6 x 10^^, 
8 X 10^^ and 2 x 10^"^ GeV^^ Fig.|4]shows the distribution of the combined probability for signal and back- 
ground events, the selection efficiency and the signal-over-background ratio as a function of the combined 
probability value, as well as the signal selection efficiency as a function of muon momentum. 

There are two main effects on the muon momentum distribution in selected events. First, the efficiency 
of the selection on the combined probability is not flat with the muon momentum. Therefore, a cut on 
this variable introduces an inefficiency at the lower edge of the distribution and a subsequent bias towards 
higher momenta, see Figure life). This inefficiency increases with the value of the probability cut applied. 
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Figure 3: Discrimination variables used in the combined likelihood function: (upper left) 6 acoUinearity, 
(upper right) sin(0Ei^iss missing energy direction, (lower left) thrust of the two muon system and 
(lower right) distribution of the variable A (see text) 

The inefficiency and the bias increase also when the momentum resolution degrades. Fig. [5] shows the 
same distributions for 5pt/p^ = 2 x 10^"* GeV^\ Fig. [Sjc) shows a bias towards higher momenta. This 
effect is accounted and corrected for in the fits performed for signal+background (see Figure [12] (b)) Both 
beamstrahlung and momentum resolution introduce a smearing of the upper momentum edge. Both effect 
have a potential impact on the statistical accuracy and the bias in extracting the SUSY particle masses from 
a fit to the reconstructed momentum distribution, as discussed below. 

The Signal-over-background ratio depends also on the beam polarisation. Fig. |6] shows the efficiency and 
S/B as a function of the probability value for different polarisation options. 

Table [U lists the number of signal (S) and background (B) events, the selection efficiency £ and the S/B 
ratio for different values of the probability cut, momentum resolution, polarisation and time stamping values. 
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Figure 4: (a) (left panel) distribution of the combined probability variable for signal events (blue) and back- 
ground events (red); 5pt/Pt^ = 2 x 10^^ GeV^\ (b) (middle panel) efficiency and S/B as a function 
of the probability value without polarisation, (c) (right panel) selection efficiency for a probability 
cut larger than 0.8 as a function of muon momentum. 
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Figure 5: Same as Fig,|4]for 5pt/p^ = 2 x 10 ^ GeV ^ In (c) the deformation of both the lower and the 
upper end of the specti^um after selection cuts is visible. 

3.3 Smuon and neutralino mass determination 

The smuon and neutralino masses are extracted from the position of the kinematic edges of the muon mo- 
mentum distribution, a technique first proposed for squarks ||9l, then extensively applied to sleptons ifTOl : 
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Eh,l = ^| 1-^ I I l±Vl-4^ I (1) 
The smuon and neutralino masses depend on the beam energy \/?/2 and the kinematic edges Eh l as: 

^2\l/2 / 1/2 



— 1 - -1— z— ^ and m^o = mr,± 1 -= (2) 
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Figure 6: efficiency and S/B as a function of the probability value for different polarisation options, (a) no 
polarisation (b) 80% e^ polarisation and (c) 80% e^ + 60% e+ polarisation 



Prob cut 


5pt/Pt^ 


pol 

{e-/e+) 


BX 


N ■ 


Nbkg 


£ 


Nsig/Nbkg 


0.80 





0/0 





1315 


2937 


0.93 


0.45 


0.80 


2 X 10-^ 


0/0 





1319 


2984 


0.93 


0.44 


0.80 


4 X 10"^ 


0/0 





1319 


2953 


0.93 


0.45 


0.80 


6 X 10"^ 


0/0 





1318 


3098 


0.93 


0.43 


0.80 


8 X 10"^ 


0/0 





1317 


3316 


0.93 


0.40 


0.80 


2 X 10-"* 


0/0 





1318 


4033 


0.93 


0.33 


0.80 


2 X 10"^ 


80/0 





1319 


1381 


0.93 


0.96 


0.80 


2 X 10"^ 


80/60 





1319 


1180 


0.93 


1.11 


0.80 


2 X 10"^ 


80/60 


5 


1317 


1271 


0.93 


1.04 


0.80 


2 X 10-^ 


80/60 


20 


1299 


1301 


0.91 


1.0 


0.90 


2 X 10"^ 


0/0 





1285 


2619 


0.91 


0.49 


0.90 


2 X 10"^ 


80/0 





1285 


1179 


0.91 


1.09 



Table 1: Scalar muon selection: number of signal, Ngig, and background, Nbkg, events for 2 ab^ of inte- 
grated luminosity, selection efficiency, e, and signal over background ratio, Ngig/Nbkg, for different 
probability cut, momentum resolution, polarisation and time stamping values. 



where Eh and Ei are the high and low momentum edges of the muon momentum distribution. This shows 
that an accurate measurement of the shape of the luminosity spectrum must be achieved and the value of 
masses extracted from the momentum spectrum are correlated. We extract the p.R and masses from a 2- 
par fit to the reconstructed momentum distribution. The fit is performed with the MiNUlT minimisation 
package ifTTTl . We model the momentum spectrum according to ([Hi, where ^/s accounts for beamstrahlung 
and ISR effects, as discussed below. Momentum resolution is included through a parametric smearing of the 
Pt distribution for the analysis performed at generator level or full tracking for simulated and reconstructed 
events. The fit also accounts for the correlations between the fiR and Xi masses. To investigate the different 
contributions to the statistical uncertainty on the smuon and neutraUno masses, several fits are performed by 
changing the input conditions. 
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3.3.1 Energy spread and ISR 



We study the contribution of the centre-of-mass energy spread to the statistical accuracy of the fit. There are 
three sources of energy spread: the momentum spread in the linac, which gives a ~7.5 GeV Gaussian smear 
on y/s for the CLIC parameters, beams trahlung, which contributes a long tail and initial state radiation (ISR); 
the first two are induced by the machine and we shall refer to them collectively as "luminosity spectrum". 
We estimate the contribution of the luminosity spectrum to the statistical accuracy on the masses and of the 
knowledge of its shape to the mass accuracy and bias. We use the luminosity spectrum obtained from the 
GuineaPig |[T2l beam simulation for the 2008 CLIC parameters. First, we compare the results of the fit 
for i) events generated without luminosity spectrum spread at ^/s = 3 TeV, ii) events in the main peak of the 
luminosity spectrum, 2950< ^/s < 3020 GeV and iii) all events with ^/s > 2500 GeV. In all these cases we 
apply a loose signal selection and we assume no resolution smearing for the muon momentum. Even without 
the luminosity spectrum contribution, the sum of the energies of the colliding electrons extends to energies 
significantly below the nominal ^/s due to QED effects. We model the ISR spectrum by an approximate 
solution to the Gribov-Lipatov equation, proposed in |[T3l . In the formula we leave free the T] parameter and 
the fraction of events off the full energy peak. We determine them by a fit to the ISR spectrum of Pythia 
signal events (see Figure |7]l. The resulting function is used to fold the ISR contribution in the shape of the 




Figure 7: Centre-of-mass energy distribution including (a) ISR and (b) ISR and beams trahlung. The points 
represent the simulation and the lines the functions used for describing their distribution in the 
mass fit. 

muon momentum spectrum used in the mass fits. Fig. [8] shows the effect of ISR and ISR -i- beamstrahlung 
on the signal muon momentum spectrum. In order to assess the effect of the knowledge of the luminosity 
spectrum on the mass measurement accuracy, we consider the luminosity spectrum obtained from CALYPSO 
for simulated signal events and we model it using the parametrisation proposed in |IT4l. This parametrisation 
has two components: a core, which we assume to be Gaussian, and a tail. We perform a to the luminosity 
spectrum with five free parameters: the width of the Gaussian core, two parameters describing the tail shape 
and two normalisation coefficients. The result of the fit is shown in Figure|7J Then, we compare the results of 
the mass fit when we use the fitted parameters of the luminosity spectrum parametrisation to those we obtain 
by varying these by ±15% of their values in a fully correlated way. This change of parameters corresponds 
to a change of the average ^/s value by ±2x10^^. The mass and statistical uncertainty of the smuon change 
by ±0.8 GeV and ±15%, respectively, and that of the neutralino by ±1.6 GeV and 10%, respectively. The 
actual accuracy on the determination of the shape of the luminosity spectrum will need to be assessed from 
a detailed study of observables such as the electr^on acoUinearity in Bhabha events 115,1 . but are expected 
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Figure 8: Signal muon momentum spectrum with no ISR/FSR nor beamstrahlung effects (Ught grey), ISR 
and FSR only (grey) and also beamstrahlung effects (black) showing the progressive smearing of 
the upper kinematic edge. 
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to be not larger than those assumed here. Fig. |9] (left) shows the fitted muon momentum distribution for 
events with 2950 GeV < ^ < 3000 GeV and (right) for events with 2500 GeV < < 3000 GeV. Results 
are summarised in Table |2l The fitted masses are in agreement with those generated mp± = 1109 GeV and 
m5>o=554 GeV, within statistical uncertainties. 
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Figure 9: Fits to the signal muon momentum spectrum for two selections on y/s. 



3.3.2 Muon photon radiation (FSR) 

A source of resolution loss is photon radiation from muons. At 3 TeV, in about 15% of the events the muon 
radiates a photon. A fit to the muon momentum distribution for signal events applying only a loose selection, 
probability cut=0.5, a cut on the centre-of-mass energy, and without momentum resolution smearing leads 
to a small increases of the uncertainty on the neuti'alino mass, but a shift on the mass value, (see Table |2l). 
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3.3.3 Event selection systematics 



The signal selection cut may introduce a bias on muon momentum distribution which propagates on the 
result of the fit to the smuon and neutralino masses. In order to study the effect of this cut, we fit the 
muon momentum distribution for signal events with a momentum resolution smearing and two different 
probability cuts, in the range 0.8 to 0.99, For a cut at 0.99 m^± = 1 127.6 ± 3.5 GeV and m^o = 557.6 ± 1 .7 
GeV. For a cut at 0.8 m,-,± = 1104.6 ± 3.0 GeV and m^o = 560.0 ± 1.6 GeV. For the events 'selected with a 
cut of 0.8 the fitted masses are in agreement with those generated, while for the tighter cut at 0.99 results 
are significantly biased. This could be eliminated by applying an efficiency correction which could carry 
systematic uncertainties. Therefore, for this analysis we adopt a selection cut at 0.8, which appears safe both 
in terms of signal-to-background ratio and signal bias. 



3.3.4 Muon momentum resolution 

Next, we estimate the contribution of the muon momentum resolution on the accuracy of the the masses 
coming from the fit. In multi-TeV collisions there is no equivalent of the Higgstrahlung e^e^ — > H^Z^ — 
Xt^i^ , {i = e, ji) process, which sets a strict requirement for momentum resolution at lower y/s values. 
Reactions such as smuon production in SUSY and ji+pL- in the SM nH can provide useful guidance 
on the track momentum resolution requirements at high energies. We express the resolution in terms of 
5pt/p^, where pt is momentum component in the plane normal to the beam axis. We perform the mass fit for 
signal events fulfilhng a loose selection and 2500 GeV < y/s < 3000 GeV assuming different momentum 
resolution values: Spt/pt^ = 0, Spt/pt^ = 2 x 10"^ GeV"', Spt/pt^ = 4 x 10"^ G&V \ Spt/pt^ = 6 x 
10-5 GeV-i, Spt/pt^ = 8 X \0-^ GeV"' and 5pt/pt^ = 2 x 10""^ GeV"'. Fig. [JO] and [H] show the fits 
to the signal muon momentum distribution for various momentum resolution values. The smuon and 
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Figure 10: Fits to the signal muon momentum spectrum for momentum smearing of (a) 5pt/pt^ = 2 x 
10-5 GeV-i and (b) 5pt/pt^ = 4 x 10-^ GeV-^ 

neutralino masses are in good agreement with the generated masses. The uncertainty on the masses starts 
being significantly impacted from the momentum resolution when 5pt/pt^ is larger than 5 x 10-^ GeV-^ 
(see Table 111). 
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Figure 11: Fits to the signal muon momentum spectrum for momentum smearing of (a) 5pt/Pt 
and (b) 5pr/p? = 2 x 10"^ GeV^^ 
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3.3.5 Background subtraction 

The cross sections for the SM processes which can lead to the same final state as the signal are one to two 
orders of magnitude larger compared to that of the jUr signal, in absence of beam polarisation. In order 
to assess the impact of the background on the statistical accuracy for the extraction of the /i« and Xi masses 
we repeat the analysis to the momentum distribution with both signal and background events. The W+IV- 
background is modelled using an "ARGUS" function ifTTl in the range >200 GeV and a first order poly- 
nomial in the range 100 GeV< <200 GeV. The other backgrounds are modelled using a polynomial 
distribution. These functions are fitted on the momentum distribution of background events passing all the 
selection cuts and used to subtract the estimated background contribution from the signal + background mo- 
mentum distribution. After background subtraction the signal distribution is corrected to take into account 
the momentum dependent selection efficiency. The fit is performed on the background-subtracted momen- 
tum spectrum. Fig. [12] shows the muon momentum distribution for signal and background events before (a) 
and after (b) background subtraction. Events are selected with a probability cut of 0.8 and the background 
is scaled assuming a 80 % electron beam polarisation. Fig. [13] shows the muon momentum distribution for 
background-subtracted events assuming (a) no polarisation and (b) both electron and positron polarisation. 
The polarisation of the electron beam only (option ii)) allows us to improve the measurement of the smuon 
and neutralino masses by 44 % and 59 % to a relative statistical accuracy of 0.8%. Adding positron beam 
polarisation (option iii)) further reduces these uncertainties to 0.6% and 0.5%, respectively (see Table [2]). 
Background rejection by the use of polarised beams is far superior compared to what can be achieved using 
tighter cuts in absence of polarisation, as shown by a comparison of the results obtained with a 0.8 probabil- 
ity cut and electron polarisation to those for a tighter cut at 0.9 for unpolaiised beams in Table|2] A dedicated 
energy scan of the smuon pair production threshold can further improve the measurements of these masses, 
also reducing their correlation. 



3.3.6 77 -^hadrons Background 

In e^e^ collisions a high rate of 77 collisions arises from photons radiated in the electro-magnetic inter- 
actions. On average there aie about 3.3 77 — )• hadrons per bunch crossing (BX). The products of the 77 
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Figure 12: (a) Muon momentum spectrum for signal + background events with highlighted the different 
components and the fitted background shape, (b) fit to the muon momentum distribution for 
background-subtracted events. Simulation assumes 80 % electron polarisation, momentum reso- 
lution 5pt/pt^ = 2x 10^^ GeV^' and selection cut value of 0.8 



interactions overlap with those from the interactions under study. At CLIC, the 312 bunches of a train, 
separated by 0.5 ns, generate a significant number of extra particles which are superimposed to the products 
of the main e^e^ events and degrade the quality of the measurement of its properties [il8]. To estimate 
the contribution of this background to the uncertainty on the smuon and neutralino masses, particles from 
77 — )• hadrons background are overlayed on signal and SM events, assuming a detector time stamping ca- 
pability con^esponding to the integration of 5 BX and 20 BX. In this analysis the main effect is the change 
in the efficiency of the signal selection. The normalised signal-to-background ratio, S/B, probabilities of 
the discriminating variables, as well as the combined probability Prob are computed for a detector resolu- 
tion: 5pt/p^ = 2 X 10^^ GeV^'. We find that for the integration of 5 BX, the selection efficiency remains 
virtually unchanged at 0.93, while for 20 BX it becomes 0.91. 

3.3.7 Full Simulation and Reconstruction 

Finally, we repeat the analysis using fully simulated and reconstructed signal events. The beamstrahlung 
effects on the the luminosity spectrum are included. The simulation is performed using the Geant-4- 
based lEl Mokka program ||20l with the CLICOl-ILD detector geometry, which is based on the ILD 
detector concept being developed for the ILC. 

Events are subsequently reconstructed using the Marlin reconstruction program [21]. Figure [14] shows 
the measured momentum resolution 5pt/pt^ obtained for muons in signal events. The masses and accuracies 
from the fit to the fully simulated and reconstructed events, (1118.4± 3.0) GeV and (569.1 ± 1.5) GeV, 
agree with those obtained at generation level with 2x 10^^ GeV^^ momentum smearing (see TableO. 



13 




Figure 13: Fit to the muon momentum distribution for background-subtracted events. Simulation assumes 
(a) no beam polarisation and (b) 80 % electron and 60 % positron polarisation, momentum reso- 
lution 5pt/pt^ = 2 X 10^^ GeV^^ and selection cut value of 0.8 



3.4 Summary 

This study allows us to draw some conclusions on the potential of a 3 TeV CLIC collider in SUSY spec- 
troscopic measurements and some of the requirements on the detector and the beams. Because of the tiny 
production cross section in the chosen high-mass scenario, background subtraction is the dominant source 
of statistical uncertainty. Electron beam polarisation at ~ 80 % gives an equivalent luminosity gain of a 
factor of six and is essential to recover precision. Positron polarisation is desirable, since it gives an addi- 
tional gain of a factor of two in equivalent luminosity and it also allows us to disentangle the contributions 
of fiL and fiR. Smuon and neutralino masses of 1108.8 GeV and 554.3 GeV, respectively can be extracted 
from the muon kinematics, in events with two oppositely charged muons and missing energy, with a relative 
statistical accuracy ~ 0.5 % with 2 ab^' of integrated luminosity and both beams polarised. In addition, the 
signal production cross section of 0.7 fb can be determined with a relative statistical uncertainty of 2.0 %. 

Since a major source of smearing of the kinematic edges of the muon momentum spectrum is beam- 
strahlung and ISR, the track momentum resolution does not appear to be critical for the measurement of the 
smuon mass, as long as a resolution 5pt/pt^ < 5 x 10^^ GeV^' can be achieved, though it remains important 
for the neutralino mass. It is important to have a good control of the luminosity spectrum and desirable to 
limit the beamstrahlung not significantly beyond that corresponding to the 2008 CLIC parameters. 

Finally, the effect of the overlay of 77 — > hadrons events from machine-induced background does not lead 
to any significant degradation of the signal selection efficiency for a detector with time stamping capability 
of 10 ns. 
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(a) 5pt/pt^ Resolution from Fully Simulated and Reconstructed Events 



1107.14± 2.81 
560.08 ± 1.49 
934.59 ,E,= 181.54 



DST:\(S>2500 

: 6 P,;P,'=2.10 '' ; polO 
Fil: entries = 2447 ,x2= 46.33 




(b) Momentum Spectrum from Fully Simulated 
and Reconstructed Events 



Figure 14: Validation using fully simulated and reconstructed events for the CLICOl-ILD detector. (Left) 
Distribution of the difference between the generated and reconstructed pt of muons normalised to 
the squared pt (5pt/p^), after full simulation and reconstruction. The width of the fitted Gaussian 
curve is 1.8 x 10^^ GeV^^ (Right) Fit to the signal muon momentum spectrum. 
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Table 2: Summary of the results of the fits to the smuon and neutralino mass for various assumptions on 
track momentum resolution, beamstrahlung, polarisation and number of bunch crossings integrated 
in one events. The results obtained on signal only (S) at generator level are also compared to those 
from full simulation and reconstruction and signal+background (S+B) fits. 







Data 


Pol 


BX 


(M±aM) (GeV) 


(xlQ-^ GeV"') 


(GeV) 


Set 


(e-/e+) 








0. 


2950 


S 


0/0 





1106.3±2.9 


558.8 ± 1.3 


0. 


2500 


S 


0/0 





1098.8± 2.6 


555.4 ± 1.2 


0. 


2500 (ISR only) 


S 


0/0 





1109.2± 3.2 


555.4 ± 1.2 


0. 


2500 


S (No FSR Cor) 


0/0 





1095.3± 3.2 


557.7 ± 1.3 


2. 


2500 


S 


0/0 





1104.6± 2.9 


560.0 ± 1.7 


2. 


2500 


S (G4+Reco) 


0/0 





1107. 1± 2.8 


560.1 ± 1.5 


4. 


2500 


S 


0/0 





1102.8±2.9 


557.2 ± 2.8 


6. 


2500 


S 


0/0 





1098.8± 3.1 


559.1 ± 3.6 


8. 


2500 


S 


0/0 





llOl.Oib 3.4 


564.2 ± 4.0 


20. 


2500 


S 


0/0 





1107.5±4.2 


575.7 ± 5.3 


2. 


2500 


S+B (0.8) 


0/0 





1107.5±15.5 


542.5 ± 11.3 


2. 


2500 


S+B (0.9) 


0/0 





1107.5±14.4 


551.2 ± 12.0 


2. 


2500 


S+B (0.8) 


80/0 





1107.7± 8.7 


542.6 ± 4.6 


2. 


2500 


S+B (0.8) 


80/60 





1118.5± 6.1 


551.3 ±3.0 


2. 


2500 


S+B (0.8) 


80/60 


5 


1105.7± 6.3 


549.4 ± 3.9 


2. 


2500 


S+B (0.8) 


80/60 


20 


1113.2± 6.8 


550.3 ± 3.4 
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